Zigzag graphene nanoribbons without inversion symmetry 
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Graphene on a substrate will suffer an inversion-symmetry-breaking (ISB) lattice potential. Tak- 
ing electron-electron interaction into account, we study in this paper the possibility of half- met allicity 
and noncollinear (NC) magnetic phase for graphene zigzag nanoribbons without inversion symme- 
try. At half-filling it is found that half-metallic(HM) state can be achieved at an intermediate 
value of the ISB potential due to its competition with the electron-electron interaction. Away from 
half-filling, the phase diagrams of doping versus ISB potential for different ribbon width are given, 
where the regimes for the HM states and NC magnetic state are clearly indicated and discussed. For 
ribbons with perfect edges, we predict a topological transition between two HM states with different 
magnetic structures, which is accompanied by an abrupt transition of electrical conductance along 
the ribbon from 2e 2 /h to e 2 /h. 

PACS numbers: 72.25.-b, 75.75.-c, 73.22.-f, 71.10.Hf 



I. INTRODUCTION 

Graphene sheet and its related nanostructures have at- 
tracted considerable attention owing to their remarkable 
electronic and structural properties [1-0] and their possi- 
ble applications in graphene-based spintronicsfj-Q- The 
free standing graphene crystallite lacks an energy gap 
and has a Dirac-cone structure due to its negligible spin- 
orbit interaction. To make graphene be used as a semi- 
conductor, several proposals have been made to open up 
a tunable gap in its electronic spectra. One is to con- 
sider epitaxial graphene on the top of a substrate which 
breaks the inversion symmetry [7H14I] . An alternative 
strategy is to construct periodic structures such as an- 
tidot lattices [lji [lg|, graphene ribbons [13, [3, graphene 
with regular patterns of hydrogen- covered regions |19ll2q| 
or with patterned defects [2ll.[23|. 

The zigzag terminated graphene nanoribbon(ZGNR) 
has attracted more attentions because it presents a band 
of zero-energy modes which is the surface states living 
near the edge of the nanoribbon [23l426| | . According to 
Stoner criterion, the ground state has a ferromagnetic 
instability due to electron-electron interactions and will 
then lead to an antiferromagnetic(AF) structure where 
each edge is ferromagnetically polarized but coupled with 
each other antiferromagneticallv|23l |27|. Most theoret- 
ical studies have focused on prope rties of the neutral 
ZGNR. Several related works [28h3(j| have considered the 
carrier density away from half-filling and the stable non- 
collinear(NC) canted magnetic states are predicted in the 
low doping regime. 

Application of a transverse homogeneous electric fields 



to ZGNR, half-metallicity[3l|can be realizedfH [Hj]. 
This means that a small longitudinal source-drain field 
could be applied to generate fully spin-polarized currents, 
i.e., electrons with one spin orientation is metallic while 
electrons with the other orientation is insulating. Some 
later works |34l - l42l | verified the promising applications for 
the future spintronics. 

In this paper, we present a theoretical study on the 
possibility of the half-metallicity and NC magnetic states 
in the ISB graphene nanoribbons. It is found that the ISB 
potential together with the carrier density affect dramat- 
ically the magnetic ground state and the band structure 
in ZGNR, leading to phase transitions between a series 
of phases including HM states. 

The paper is organized as follows. In Sec. II we give 
the model Hamiltonian and its mean- field treatment. In 
Sec. Ill we discuss the magnetic structures of the ground 
states. In Sec. IV we investigate and discuss the HM 
states and correspondingly give the band structures as 
well as the phase diagram. In Sec. V, we summarize our 
results. 



II. MODEL 

We assume that the bulk graphene is subjected to a 
staggered sublattice potential, which breaks the inversion 
symmetry of graphene. This is the general situation of 
the graphene on a substrate. The model Hamiltonian can 
be written as follows, 
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Here t is the nearest-neighbor hopping integral, and 
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FIG. 1: (Color online) Representative magnetic structures 
for an ISB ZGNR with width W = 10. (a) AF ground 
state for V — 0.10 and Sn — 0.1, where m L = 0.091, 
m R = 0.129.(b) FM ground state for V = 0.04 and Sn = 0.2, 
where m h = 0.066, m R = 0.135. (c) Fb ground state for 
V = 0.10 and Sn = 0.4, where m L = 0.0009, m R = 0.125. 
(d) NC ground state for V = 0.02 and Sn = 0.046, where 
m L = 0.129, m R = 0.133 and 9 — 68°. Here m L , m R are 
the magnitudes of the spin polarizations on the left and right 
edges, respectively, whereas 9 is the relative orientation angle 
between them. The circles(dots) denote the triangular sublat- 
tice A(B). The dashed rectangle denotes the unit cell, which 
is periodically repeated along the vertical direction. 



FIG. 2: (Color online) (a): The relative orientation angle 9 
between the spin directions on the two edges in the NC canted 
states as a function of doping at different V for an ISB ZGNR 
with width W = 10. The metastable NC states are denoted 
by the cross marked points, (b)-(d): Energy comparison of 
the three magnetic states(NC, AF and FM) as a function of 
doping for the same ribbon at different V: (b) V — 0, (c) 
V = 0.02, and (d) V = 0.03. The energies are all shifted by 
that of the AF state E AF at the corresponding doping value, 
so as a result the triangles at the zero energy denote the AF 
state. The filled circles(squares) denote the energy difference 
between the the NC(FM) and AF states. 



U the on-site Coulomb repulsion energy. Vi is the stag- 
gered sublattice potential with Vi = V(— V) on sublat- 
tice A(B), and n is the chemical potential. Although 
the substrate potential breaks the inversion symmetry, 
the Hamiltonian is still invariant under PI operation, 
where P is the particle-hole transformation^ ci CT — ^Cj^, 
with r] = 1(— 1) if i belongs to A(B) sublattice) operator 
and I is the space inversion operator, with the inversion 
center chosen as the center of one of the central hexagons 
of the ribbon. In the following, energy is measured in unit 
of t^2.8eV. 

To take into account the possibility of the noncollinear 
spin polarization, the Hubbard term is so decoupled that 
the mean-field Hamiltonian can be written as, 
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where the electron spin polarization vector and the 
charge density are given by m^ = (l/2)(cj<rci), and 
rii = (ctcj), respectively, with cr = (<T x ,a y ,a z ) the Pauli 
matrices and c\ — (cL,cL). In the whole paper, we 
choose U/t = 1.0, which is an appropriate value con- 
sistent with the first-principle calculations J43|. For this 
intermediate U , the decoupling process introduced is be- 



lieved to be reliable and effective. Starting from an initial 
random spin and density configuration we determine the 
above parameters for a sample with a strip geometry self- 
consistently by the standard iteration method. 



III. THE MAGNETIC STRUCTURES 

Similar to ZGNR with the inversion symmetry, the 
spins at the edges of the sample will be polarized and 
the solutions can be classified into several categories ac- 
cording to the coupling between the spin polarizations at 
the two edges. For different ISB potentials and differ- 
ent carrier densities, one can find the following ground 
state solutions: the AF state and ferromagnetic(FM) 
state, where the spin orientations at the two edges are 
antiparallel and parallel to each other respectively; the 
so called Fb state|28|], where only one edge is spin polar- 
ized; and the non-collinear(NC) canted state where there 
is a relative orientation angle between the two spin 
directions. The corresponding magnetic structures are 
schematically shown in Fig. [TJ where the effective doping 
value Sn = (1 - (N)) x 2W with W the width of ribbon 
and (N) the average electron number per site, is mea- 
sured as the density departure from half-filling per unit 
cell of ribbon. 

The NC canted state is found to be the ground state 
only in the very low doping level and small V. With the 
increase of V or ribbon width W, the doping regime in 



which NC is favorable will become narrower and finally 
disappear (see the phase diagram in the next section). In 
Fig. [2a) we plot as a function of doping at different 
V. For a very small ISB potential, 9 gradually decreases 
from 180° to 0°, whereas for a finite ISB potential, NC 
solutions can be found only above a critical doping level 
which also increases with V. Some NC solutions are actu- 
ally metastable and have a higher energy than the AF so- 
lutions. In Fig. [^b)- (d) we compare the energy of the NC 
state with that of the AF and FM states, which clearly 
show the stability of the NC ground state at low-doping 
levels. 

These NC canted states are found to be always 
metallic. The representative band structure is shown 
in Fig. EJa), which is characterized by four counter- 
propagating current-carrying states at the chemical po- 
tential, leading to a quantized electric conductance 2e 2 /h 
along the ribbon with perfect edges. These four partially 
polarized current-carrying states (see Fig. [3jc)-(d)) are 
not edge ones but extended along the transverse direc- 
tion even for a wider ribbon(see Fig.[3jb)), which can be 
due to the relatively strong coupling between two edges. 
Detailed calculation shows that the NC canted state for 
V = case shares the same picture, but it has a spe- 
cial property which is not preserved by V ^ case that 
the averaged spin carried by the four propagating states 
is exactly parallel or anti-parallel to the vector sum of 
the two edge polarizations m R + m L (see the caption of 
Fig.©. 



IV. HALF-METALLIC STATES AND THE 
PHASE DIAGRAM 

A. Half-filling 

1. Classification of the ISB graphene nanoribbons 

There are two competing processes in the ISB graphene 
nanoribbons. One is electron-electron interactions, the 
other is the ISB external field. In a neutral ZGNR with 
the inversion symmetry (V = 0), the ground state has 
the AF magnetic structure, which is a consequence of the 
flat subbands at zero energy. A large magnetic moment 
emerges around the edge sites even for an infinitesimally 
small Coulomb U, inducing a finite gap [44], [45( as a re- 
sult. On the other hand, if we neglect electron-electron 
interaction, application of an ISB external field will mix 
the wave functions at the two sublattices and will then 
open up a bulk band gap in neutral ZGNR. Since either 
of the processes prefers to open up a gap, it is quite un- 
usual that when both of them are present in a graphene 
ribbon, a state with a negligible gap can be achieved. 

For an ISB graphene ribbon, it is found that there 
exists a ribbon- width dependent threshold value V c of V, 
beyond which the system becomes paramagnetic(PM). 
When V < V c , the ground state at half- filling (/x = 0) 
has an AF magnetic structure which preserves the PI 




0.10 


. (b) 


— mM 2 


- 


0.08 


I, 


— mM 2 


- 


0.06 


It 


- Ms)f 


- 


0.04 


i 


ffCyf 




0.02 
ono 


»«! Ill, ,ii 





0.02- 


(c) 






0.00- 


ir 






0.02- 


-■- s z< k ,) 

-•- s,(ig 






i 


- SJk,) 




0.04- 




— S ,(K) 






20 40 


60 


30 



0.02 
0.00 


(d) 


0.02 
0.04 


fMlltW 1 "' 1 


— s ,( k ,> 

-•- S x( k 2> 

- S x( k ,) 

— S x( k «) 


- 




) 20 40 


60 


80 



(e) V=0.01 




(f) V =0 s=s 3 






s„=s. 



FIG. 3: (Color online) (a) Typical band structure for a NC 
metallic state at V = 0.01 and on = 0.024 for ZGNR with 
width W = 40, where the angle — 80° and m R ~ m L = 
0.138. There are four counter-propagating current-carrying 
states ki(i = 1,2,3,4) at the Fermi level. Correspondingly, 
the electron density and spin density distributions for the 
four states are shown in (b)-(d). If the spin axis is so cho- 
sen that the edge polarizations can be written as ni R(L) = 
(±sin(0/2), 0, cos(0/2)), the averaged spin carried by the four 
states can be expressed as Si = mi(cos(9i),0, sin(6i)), where 
mi = m 4 = 0.400, m 2 = m 3 = 0.353, and 0i = 04 = -72.3°, 
02 — 03 — 51.5°. As a comparison, the result for the NC 
canted state with the same doping Sn — 0.024 for V — 
case is given as follows: 8 — 64°, m B — m L = 0.138, 
mi = m 2 = 0.424, rn 3 = m 4 = 0.299, and 8 X = 04 = -90°, 
02 = 0» = 90°. These vector quantities are schematically 
shown in (e) and (f). 



symmetry: 



(3) 



implying antisymmetric and equal spin polarization on 
both edges, with n R „ , the electron density at the right- 
most(leftmost) boundary sites. Therefore, the electronic 
spectra for the spin-up and spin-down components can 
always be decoupled and both has the symmetry about 
zero energy. In Fig.[4l band structure for a neutral ZGNR 
with width W = 20 for both spin directions are shown. 
In each sub-figure, the center four bands represent that 
for the edge states. Among the four bands, near k y = n, 
the red and black ones below(above) zero energy corre- 
spond to spin-up edge states on the right(left) boundary 
and spin-down ones on the left(right) boundary respec- 
tively. Upon increasing V, though the bulk gap remains 
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FIG. 4: (Color online) (a)-(e):Band structures of ZGNR with 
width W = 20 at half-filling for different ISB potential V: 
(a) V = 0.03; (b) V = 0.05; (c) V = 0.11; (d) V = 0.13; 
(e) V — 0.15. From (a) to (d), the ground state has the AF 
magnetic structure and the red and black lines denote the 
spin-up and spin-down states respectively. The ground state 
in (e) is a non-magnetic one. The denoted squares are the 
states at k y — tt which are completely localized at the zigzag 
boundary sites. Correspondingly, the spin polarizations at 
the leftmost and rightmost boundary sites are shown in (f) as 
functions of V. 



open, the gap between the two spin-down edge states is 
becoming smaller. When 0.116 < V < 0.136 = V c , the 
gap is vanishingly small(less than 10 _4 i) while spin-up 
edge states are still gapful. Thus, a HM state is realized 
in this regime. We note that the situation is quite similar 
to that studied in the HM state in ZGNR when apply- 
ing a transverse electric field 32], which also breaks the 
inversion symmetry of graphene ribbon. 

The reason of the gap opening and closing can be at- 
tributed to the interaction between the left and right 
spin-down edge states, which is similar to the gap- 
opening in finite quantum spin Hall systems 46]. Ac- 
cording to the perturbation theory, the gap due to their 
interaction can be expressed as A/t = |^*(fc*,x — 
i) , 5 L (A;*,x)|. Here the "crossing point" k* is the wavevec- 
tor where gap is opened up, and $ L (fc*, x), ^ R (k*, x) are 
the corresponding spin-down edge-state wave functions 
for the half-infinite ZGNR with left and right zigzag edges 
respectively, with x taken to be the leftmost site of the 
ribbon unit cell and i = (l,0)a. Actually, the states at 
k y = 7T are completely localized at the zigzag boundary 
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FIG. 5: (Color online)Band structure of the ground states as 
a function of edge potential for a ZGNR with width W — 20 
and V — 0.05. Here the right-most site potential V R is fixed 
to be zero, whereas the left-most site potential V^ is tuned 
from being negative to positive. 



sites. Their eigenenergies can be given analytically by, 
K„ (k v =n)=V + U(n L - l)/2 - am L (4) 



E H „(k v =Tr) = -V + U(n H -l)/2 



(5) 



satisfying E^ (k y 



-E R<r (k v ). So the edge states near 



k y = ir are well localized and their decaying lengthes 
are very small, leading to a vanishingly small overlap 
and then a vanishingly small gap. On the contrary, the 
edge states away from k y = tt are much more spread and 
have a relatively large overlap so a finite gap is opened 
up. When V > V c , the AF magnetic structure is un- 
stable and the ZGNR will undergo a phase transition 
from the AF HM state to a charge-density-wave insula- 
tor PM state which has a charge gap of magnitude of 
2V - U(n R - n L )/2(see Fig. EJe)). The above discus- 
sion is based on the result of a graphene ribbon with a 
fixed width W — 20, but detailed calculations show that 
although the threshold V c and the band gap may vary 
quantitatively a little with the width, the fact of the ex- 
istence of HM regime is not changed qualitatively. 

Generally, graphene on a substrate will inevitably 
suffer a Rashba spin-orbit (RSO) interaction H R = 
i\ J2 c \{ct x dy) z Cj, where A is the coupling strength 

<i,3> 

and dy is the unit vector along ij direction. The most 
significant effect of RSO interaction on magnetism is that 
the spin polarization will be locked to an easy plane which 
is normal to the graphene plane and zigzag line. For a 
small A, the relative angle between the two edge polariza- 
tions will be slightly deviated from 180°, and the devia- 
tion angle 69 is nearly linearly dependent of A with 56 / 5\ 
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FIG. 6: (Color online)The phase diagrams: doping value 5n 
versus the ISB potential V for a ZGNR with W = 10 for (a) 
and W — 40 for (b). The dashed line shows the sequence of 
phase transitions with increasing doping for a particular fixed 
V. 
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FIG. 7: (Color online)Two HM states for an ISB ZGNR with 
width W = 40 and V = 0.08, indicated by the two star 
marked points in Fig. [6] The left column: the HM state with 
the AF magnetic structure, where 8n — 0.1. The right col- 
umn: the HM state with the Fb structure, where 8n = 0.328. 
Accordingly, the band structures are shown in the top panels, 
and the density distributions of the edge states at the chem- 
ical potentials are shown in the middle panels, whereas the 
spin density for the ground states are shown in the bottom 
panels, respectively. 



estimated as ~ 0.5° /meV by our numerical calculation. 
In actual situation, this interaction is very small and 
can be neglected in most of applications. Furthermore, 
though RSO interaction breaks the inversion symmetry, 
it can be easily checked that this term still preserves the 
PI symmetry of our model Hamiltonian. Therefore, qual- 
itatively RSO interaction should only cause very little 
deviation from our conclusion on half-metallicity and we 
then neglect its effect in the following discussion. 



2. Manipulating graphene properties by edge potentials 

Since the HM state only exists within a limited small 
region of V in neutral ZGNR, it is meaningful if one 
can find some methods to realize half-metallicity in other 
restricted region. One way to do this is by modifying 
the edge potentials, since the edge states are localized 
around the sample boundary, and so they are sensitive 
to the variation of on-site potentials on the boundary. In 
a free-electron ZGNR system with broken bulk inversion 
symmetry, this problem with varying one edge potential 
has been studied recently in terms of valley Hall effect 
and it is found that the edge bands can be continuously 



changed by tuning the on-site edge potentials [4 7]. 

Here we examine how to manipulate ZGNR to achieve 
HM state in the presence of electron-electron interac- 
tions. Depending on the values of V L and V R , half-filling 
ZGNR can have various magnetic structures, which af- 
fect dramatically the edge bands. In Fig. [5j for a fixed 
V we give the band structure of the ground states on a 
W = 20 ZGNR as V^ is tuned from negative to posi- 
tive value with V R = left unchanged. With variation 
of V L , the system undergoes a series of phase transitions 
due to the magnetic transitions at the edges. When V L 
deviates from zero a little, the ground state still has an 
AF magnetic structure, but asymmetric(Fig.[5jd)(e)(f)). 
In a narrow regime around V L ~ 0.15, one can achieve 
an AF HM state(Fig. [5jf)). In this state, the two spin- 
down edge bands are mixed due to their interactions. 
Both of the spin-down edge bands are partially filled with 
the spin-up edge bands still gapful, leading to the half- 
metallicity of the state. For a perfect ribbon, this HM 
state is characterized by a quantized spin-down conduc- 
tance 3e 2 //i, since there exist one hole pocket and two 
electron pocket at the Fermi level and thus three pair of 
propagating modes for the spin-down edge bands (see the 
inset of Fig. iff)). 



B. Away from half- filling 

When the ISB ZGNR is doping away from half-filling, 
the HM state with the AF magnetic structure can still 
survive, except breaking the PI symmetry spontaneously. 
This is consistent with the case of doped ZGNR under a 
transverse electric field, where an AF half-metallic state 
is found to exist[30j. To have a concrete picture, we give 
the phase diagram for ZGNRs in Fig. H>] Different from 
the case at half-filling, another type of HM state with the 
Fb magnetic structure emerges at finite doping and rel- 
atively large ISB potential V. This state is HM because 
both the spin-down band are partially filled, similar to 
that induced by edge potential discussed above. When 
the ISB potential is within this regime, the system at 
half-filling is actually a PM insulator, which means that 
any small doping away from half-filling will induce a mag- 
netic transition and then lead to the phase transition 
from PM insulator to Fb half-metal. Different to the AF 
HM state which has four edge states propagating along 
the two zigzag boundaries, this Fb HM state has only 
two edge states which are counter-propagating along the 
polarized right edge(see Fig. [7]). Therefore, this Fb HM 
state is characterized by a quantized electrical conduc- 
tance with value e 2 /h, not 2e 2 /h. The phase diagram 
for a wider ZGNR(Fig. Ob)) is qualitatively the same 
as that for a narrower one(Fig. HJa)), except that there 
emerges an additional AF metallic state in the former 
case. When V = 0, our result is consistent with the 
previous work|28(. We note that in epitaxial graphene 
a molecular-doping induced metal-insulator transition is 
observed in ARPES experiment 48] . Moreover, a recent 
first-principle calculation on epitaxial graphene on SiC 
has found the evidence of HM state [49|. 

With the increase of Sn, one can realize in order the 
following phases from AF insulator: AF half-metal, FM 
metal, Fb half-metal, Fb metal and PM metal(see the 



dashed line in Fig. EtR)). Since graphene on different 
substrates suffer different ISB potential |7HI4|. upon dop- 
ing graphene with fixed ISB potential will undergoe phase 
transitions between HM states and normal-metal states. 
Particularly, graphene upon doping can change by a first- 
order transition from one HM state to another, since 
there exists a phase boundary between the two types of 
the HM states. This process may be served as a kind of 
channel switch, which "close" the propagating channel 
on the left edge by spin depolarization upon doping(see 
Fig. EJ. This is expected to have some important appli- 
cations in spin related transport. 



V. SUMMERY 

In summery, the magnetic and band structures of an 
ISB graphene ZGNR have been investigated. At half 
filling, half-metallicity is found to be realizable for an in- 
termediate ISB potential because of its competition with 
the electron-electron interactions. For doping away from 
half-filling, phase diagrams for varying ribbon width are 
given and another type of HM state is achieved with dif- 
ferent magnetic structure. Spin canted states are found 
to exist only at very low doping levels and low ISB po- 
tential. Due to the importance of half-metallicity in spin 
related transport, these results may have important ap- 
plications in graphene based spintronics. 
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